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Abstract 

The elementary excitations in superfluid liquid ^He named rotons 
have an unusual dispersion curve. The energy is an approximately 
quadratic function of (p — po), the difference between the magnitude 
of the momentum p and a characteristic value pq. As a result, while 
for p > Pq a roton has its (group) velocity parallel to its momentum, 
when p < Po the velocity and momentum are antiparallel. When 
p = Pq, the roton has non-zero momentum but zero velocity. These 
kinematic properties lead to unusual trajectories when rotons scatter 
or experience external forces. This paper examines this behavior in 
the classical (ray optics) limit, where the roton wavelength is small 
compared with all other dimensions. Several experiments illustrate 
these effects. The examples are interesting in themselves, and also 
offer unconventional pedagogical possibilities. 



1 Introduction 



1.1 Background 

The concept of elementary excitations in superfluid liquid helium (^He) was 
introduced by L. D. Landau in 1941. [jl| The essential idea of Landau is that 
the low-lying excited states of the macroscopic system have energies E {Eq 
is the ground state energy.) and momentum P which can be expressed in 
the form 

* Author's address: H. M. Randall Laboratory, University of Michigan, Ann Arbor, MI 
48109-1120, USA. 
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E = Eo + Y.^{p)e{p) 

p 

P = ^n{p)p. 

p 

That is, the system of strongly-interacting particles (He atoms) is replaced 
by a collection of non (or weakly) interacting "elementary excitations" or 
"quasiparticles." In Landau's paper of 1941, he suggested a form for the 
dispersion curve (the relation between energy e and momentum p) for these 
entities. For an isotropic system, such as a liquid, the energy can depend 
only on p, the magnitude of the momentum. The 1941 dispersion curve has 
two branches: There is an "acoustic" branch with e{p) = cp, representing 
longitudinal phonons (c is the speed of sound.), and a second branch with 
s{p) = + A. Landau called these elementary excitations "rotons."[| 

The situation seems reminiscent of that in a solid, with its acoustical and 
optical branches. There is, however, an important difference: The liquid 
has complete translational invariance (not just against finite translations). 
Thus the momentum of the elementary excitations is not just a quantum 
number like crystal momentum, but is the usual physical momentum. Landau 
evaluated the two free parameters A and fi by fitting measured values of the 
specific heat of the hquid, and used the model to calculate the normal and 
superfluid fractions as functions of the temperature. 

In his brief 1947 paper, |^ Landau^ modified the relationship e{p) in 
order to account for some experimental data on the propagation of second 
sound. The curve he drew in 1947 agrees remarkably well with modern 
determinations of e{p), using inelastic neutron scattering. The two branches 
are replaced by a single branch which has a linear part, a maximum, and 
a minimum: We call the elementary excitations near the origin "phonons" , 
those near the minimum "rotons." These two parts of the spectrum are 
dominant in determining the low-temperature properties of liquid helium. 
Donnelly has written a recent introduction to the subject with interesting 
historical perspective and some discussion of the microscopic nature of the 
roton, a subject beyond the scope of the present article. 

^Landau [Q writes that "This name was suggested by I. E. Tamm" 

^Landau's papers of 1941 and 1947 are reprinted in several books, such as |3| and [Q 
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Landau emphasized the importance of the finite slope of the dispersion 
curve at the origin — in distinction to the behavior of the Ideal Bose Gas. A 
consequence of this finite slope is a "critical velocity" , below which the fiow 
of the helium is dissipationless. In the simplest case the primary dissipation 
process would be roton emission. With rare exceptions ^ other dissipation 
mechanisms occur at velocities below that required for roton emission. 

The reasons why there are no other states near the ground state of an 
interacting Bose system are considered in the papers of Feynman in the 
1950s. 

1.2 Current interest 

Since the development of the dilution refrigerator, it has become practical 
to perform experiments on "ballistic rotons." That is, at sufficiently low 
temperature the mean free path of a roton becomes as long as the size of 
the apparatus. Rotons (and phonons) are produced by a "heat pulse" in the 
liquid or, in some cases, by helium atoms incident on the surface of the liquid 
from the (low-density) vapor above it. We will consider such experiments in 
sections | and ||. 

2 The dispersion curve 

The dispersion curve for elementary excitations in superfiuid ^He, determined 
by inelastic neutron scattering §]. is shown in figure |l]. The region near the 
minimum is usually described by the equation 

£(p) = (p-po)V2/i + A. (1) 

At low temperature and pressure, approximate values for the parameters 
are: go = Po/^ = 1-92A , A/ks = 8.62 K, = m^/Q. The function 
6{p) can be considered a Hamiltonian for the free roton, since de/dpi = Vi 
is a component of the (group) velocity of a roton wave packet. Using the 
expression (|I]), we have v = p{p — po) / /i. 
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Figure 1: Dispersion curve for elementary excitations in liquid helium, cal- 
culated from the polynomial fit in reference 0] . The roton region is near the 
minimum. 
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2.1 Statistical Thermodynamics 

We recapitulate some results which we will later find useful. Since we will 
be interested in temperatures less than or of order 1 K, we will always have 
s{p) > A >> ksT, whence exp{—e{p)/kBT) = exp(— << 1 will 
always be valid. (The symbol /3 stands, as is usual, for l/ksT.) This permits 
considerable simplification of various equations. Furthermore, all integrand 
functions of p will have contributions peaked near p = Pq. Thus, we transform 
integrals over p over the range from zero to infinity to integrals over (p —po) 
from — oo to oo. 

The expected number of rotons in a given mode p: 

{nip)) = (exp(/?£(p)) - 1)-^ exp(-/?£(p)). 
The number density of rotons: 

POO 

Up = Np/V = {2Trhy^ / exp(-/3£(p))4Vdp 

J 

The Helmholtz Free Energy: 

POO 

Fp = -kBTV/{27rhf / exp{-l3e{p))A7Tp^dp 

J 

The average "kinetic" energy per roton: 

{{p-p^f/2p) = (//^V2) 

= -d/dl3 (log exp (-I3{p - pof/2ii) p^dp^ 

The pressure due to the roton gas: 

Pp^- (dF/dV)^ = -F/V = npksT = n^/x (v^) (2) 

3 Free rotons 

It is worth repeating that the expression for the velocity has the curious 
property that v is parallel to p when p > po, but antiparallel when p is below 
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the minimum. We will refer to rotons with momentum above the minimum 
as or ordinary, those with momentum below as R~ or extraordinary. The 
surface e{p) = A is a sphere of radius po in momentum space. For an energy 
greater than A (but below the maximum) the surface has two concentric 
spherical sheets, one with p > po, on which the rotons are ordinary, the 
other with p < pq, describing extraordinary rotons. Figure |^ shows a two- 
dimensional representation. We will often denote the two "conjugate" values 
of momentum with equal energy as p and p. 




Figure 2: Surfaces of constant energy in momentum space for e = A (dotted), 
and a higher energy (solid). 



4 Rotons in a velocity field 

A somewhat elusive argument, due to Landau, yields the result that to pro- 
duce an elementary excitation of momentum p in a fluid moving at velocity 
V requires energy 

e(^^{p) = eip)+p-v, (3) 
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where e{p) is the energy of the excitation in the stationary fluid. In a slowly- 
varying velocity field, we add a term to the hamiltonian 



5e{p) = p- v{r). 



(4) 



A quasiparticle of a particular momentum magnitude has higher energy 
if its momentum points "downstream" than if it is directed "upstream". 
Thus, in thermal equilibrium, the excitations will reduce the momentum of a 
moving superfiuid. Landau uses this argument to calculate the normal-fluid 
density. In section |^, we will apply it to the problem of roton motion near a 
quantized vortex. 

We present a rather pedestrian but (we hope) straightforward derivation 
of this familiar result in Appendix |A.1| 

5 Rotons in a force field 

The roton parameters A, po, and /i can vary with position, since they are 
known to depend on the density of the liquid. The "energy gap" A, in 
particular, decreases from 8.62 K at the vapor pressure to approximately 
7.2 K at a pressure of 25 atm p|. If the roton parameters vary sufficiently 
slowly compared with the roton wavelength, the roton will follow a classical 
trajectory. Since go = Po/^ ~ 2A , we have A ~ SA; a spatial variation of A 
on a scale large compared with SA will appear like a classical potential energy 
in e{p) = Ti.{p, r). In the case in which only variation in A is considered, we 
have a Hamiltonian consisting of a "kinetic" energy T = (p — po)^/2/i, and a 
"potential" energy V{r) = A(r). 

5.1 A constant force 

An interesting illustration of the consequences of the unusual dispersion rela- 
tion is provided by considering the motion of rotons in a uniform force field, 
i.e. when A is a linear function of a single coordinate, say y. If we write 
/ = —dA/dy, then Hamilton's equations take the form 



Px 



dn/dx 







Py 



dn/dy 
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= / 

X = dH/dpx 

= (Px/p) X dT/dp 

y = OH/dpy 

= (py/p) X dT/dp. 

The equations for p can be integrated immediately when the force / is 
constant. We obtain p^, = p^i, py = Pyi + ft. On a momentum-space diagram 
such as that of figure Q the motion is parallel to the Py axis, in a direction 
given by the sign of /. For / < 0, we show the various possible cases in 
figure 1^. Incident rotons with energy e and p > po (ordinary rotons) exhibit 
three different types of trajectory, depending on the value of the angle of 
incidence. Extraordinary rotons (with p < pq) can only become ordinary. 




Figure 3: Momentum space diagram showing initial and final momenta cor- 
responding to reflection of rotons by a force in the —y direction. 

We have dT/ dp = {p — po) / /i 

Px(t) = px{0) 
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Pyit) = Py{0) + ft 

We show in figure ^ coordinate space trajectories derived from tliese re- 
sults. Tliere is a constant repulsive force in the region ?/ > 0, no force for 
?/ < 0. A cusp appears in the trajectory when p becomes equal to po? at 
which point the velocity is zero, but the momentum vector is not. The unit 
of distance has been taken as the initial "kinetic" energy divided by the mag- 
nitude of the force, so that the roton stops after penetrating unit distance 
into the region. Three rotons with momentum l.lpo have been started at 
the point (0, -1) with tangential (x) momenta 1.05, 0.95, and 0.85 times Pq. 
The trajectory of an R~ roton is similar to the last of these, but traversed 
in the opposite direction. The trajectories of rotons with different values of 
the incident momentum would be similar to these, differing only in a scale 
factor in the x direction. 




Figure 4: Coordinate space trajectories corresponding to momenta in the 
three regions shown in the previous figure. Notice the different scales in the 
two directions. 
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Such cusped trajectories were apparently first noted by Goodman 
in analyzing roton trajectories near a vortex line, a problem which we will 
discuss in section |^. Note, however, that these "snap back" trajectories (an 
apt term used in ref. [|10i) have momentum transfer 5p of order 



5p ~ P — P 
= 2(p-po), 

a quantity which is small compared with 2po. This reduces significantly the 
pressure exerted on a wall by a roton gas. This point has been emphasized 



by Maris and Cline [O . 



5.2 Pressure 

The results of the previous section imply that the pressure exerted by a roton 
gas must be much below what one would naively estimate from the fiux and 



a momentum transfer of order pq. As emphasized by Maris and Cline |Tl 
the momentum transfer must be low to agree with the exact result, given in 
equation |[ 



p, = -{dFp/dV)T = UpkBT. 

In the semiclassical case, considered here, the R~ rotons always change 
mode, never reflecting specularly; the rotons do the same, except for 
nearly glancing angles. Summing the contributions to the pressure from 
these various processes does yield the correct result. 

In quantum mechanics "Whatever is not forbidden is allowed." In partic- 
ular, specular reflection is always consistent with the conservation laws, even 
when it does not occur classically. When we consider the branching to the 
other possible flnal state, in quantum mechanics, it can not change the result. 
The pressure is independent of the details of the processes at the walls. We 
can see explicitly that there will be cancelation between specular reflection 
of rotons, which has large positive momentum transfer, and that of R~ 
rotons, which produces a similar negative contribution to the pressure. 
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6 Rotons at a free surface 



In an early paper which considered balhstic rotons, P. W. Anderson |T2 
pointed out that, since the minimum roton energy (8.6 K) is greater than the 
binding energy of an atom in the hquid {Eb = 7.2 K), a roton incident on 
the surface might eject an atom into the vapor. 

If the process is elastic, so that all the energy is transferred between the 
elementary excitation and the atom, and the surface is smooth, so that the 
component of momentum parallel to the surface is conserved, then we can 
define an "index of refraction" in the usual way. (Note, however, that an R~ 
roton and an atom would have oppositely directed tangential velocity.) 

h = Pa/p = sinO/ sin0, 

where p and 6 are the momentum and angle of incidence of the elementary 
excitation and pa and the corresponding quantities for the helium atom. 
In this section, in order to connect with experiments described in section 
we consider the entire excitation spectrum — not just the roton region. 

The index n is a single-valued function of p, but a multiple-valued func- 
tion of Pa, since an incident atom might produce as many as three different 
excitations with a given energy. Written in terms of p (or q = p/h), we have 



fi = J{2mJp^){eip)-EB] 



Figure ^ plots the refractive index versus momentum, and also shows the 
energy of the elementary excitation and a free atom of equal momentum. 
The index equals 1 where the two energies cross. 

7 Roton-vortex interaction 
7.1 Quantized vortices 

The spatially-varying superfiuid velocity field near a straight vortex line is: 

'Vsir) = K X r/27rr^ = {K/27Tr)(p. (5) 

Here r is the two-dimensional (polar) radius vector, the polar angle, and 
K is along the vortex line {z direction). The definition of the circulation k is 

K = (p Vs ■ dl 
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Figure 5: Energy of an elementary excitation, and an atom outside the liquid, 
and the "Refractive index" plotted against momentum/ h. 
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and, in the case of superfluid liquid ^He, is quantized in units of /i/m4. These 
equations are vahd outside a core region, which is of atomic dimensions. 



7.2 Quasiparticle-vortex interaction 



Goodman ||T0[ and Samuels and Donnelly [0 calculated the "mutual fric- 



tion" between normal and superfluid using the roton- vortex interaction, which 
causes momentum transfer between the rotons (a constituent of the normal 
fluid) and the superfluid vortices. 

If a quasiparticle is not too close to the vortex core, the velocity field will 
not vary very rapidly, and we can combine equations § and ^ to write: 

6e{r) = p- K, X r/27rr^, 

which can be rewritten in the form 

6e{r) = R ■ r X f) 1 
= K,Lz/27rr'^, 

if we take the z axis along the vortex line. In this form, it is apparent that 
Lz, the z component of angular momentum will be a constant of the motion, 
permitting analysis of the trajectory to be reduced to quadrature. |]14| 



We show trajectories of i?"*" and R~ rotons, calculated by numerically 
integrating the equations of motion, in figure |^. An roton with incident 
X momentum {pi)x = — 1-lpo is started from the right, an R~ roton with 
X momentum {pi)x = 0.9po from the left for various values of the impact 
parameter. Distances are measured in units of a "critical" distance 



at which the roton-vortex interaction energy is equal to the initial roton 
kinetic energy, about 16A for the rotons illustrated. 

Both the R'^ roton traveling to the left and the R^ roton traveling right 
have their momenta directed to the right. Thus, both experience a repulsive 
interaction above the (counter-clockwise circulating) vortex, where p ■ vl is 
positive. The mode-changing trajectories all have impact parameters on the 
same side of the vortex. 
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Figure 6: Calculated roton trajectories near a quantized vortex line, located 
at the origin. The curves are calculated for an incident roton with mo- 
mentum l.lpo incident from the right (solid lines), and an R~ roton with 
momentum 0.9po incident from the left (dashed lines). 
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8 The experiments 



Several experimenters performed pioneering experiments involving rotons at 
temperatures low enough so that the roton mean free path is long. 

At Ohio State University, evidence for roton production was sought in a 
study of helium atom reflection from the free surface. 



In Paris, the inverse experiment was performed. [|T6| Ballistic rotons, pro- 
duced at a pulsed heater 'evaporated' atoms from the free surface. 

These experiments have been extended and elaborated in a series of ex- 
periments at University of Exeter. In the Exeter experiments, rotons (and 
phonons) are produced at a pulsed heater and detected by a superconduct- 
ing bolometer. The experiments are performed at temperatures sufficiently 
low that the rotons may considered to be collision-free (after they leave the 
vicinity of the source), and to follow "ballistic" trajectories. We cite only a 
few of the articles here, and refer the reader to these articles for additional 



references. A rather comprehensive early article is. [|T7 



The experiments include detection and measurement of roton-roton scat- 
tering, produced by crossing two beams. 



A puzzling aspect of these experiments is the observation that only i?"*" 
rotons are produced at the heater, although there was evidence that R~ 
rotons could be produced at surfaces. 

8.1 Direct detection of R rotons 

A compelling demonstration of the odd kinematic properties of the roton 
would be evidence of oppositely-directed momentum and velocity. 

Tucker and Wyatt have recently realized a source of R~ rotons, using 



two facing heaters, and utilizing scattering of i?"*" rotons to produce R~ ones. 
They have directed these rotons at the free surface of the liquid, producing 
evaporation of single atoms. As noted in section ^ the R~ roton and the 
atom will have oppositely directed tangential velocity. 

This type of confirmation of the anomalous connection between velocity 



and momentum of the R~ rotons pO| is accomplished in the experiment 
pictured in Fig. ^ 

The evaporated He atom is detected in the "backward" direction, al- 
though the momentum component parallel to the surface is conserved in the 
transition from roton to evaporated atom. 
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A Appendix 



A.l Excitations in a moving fluid 

We present a rather pedestrian but (we hope) straightforward derivation of 
this famihar resuh. If we first think classically of viewing a system from a 
reference system in which it is moving at velocity v, each particle will have mv 
added to its momentum, which will change the momentum and kinetic energy, 
but for a system which is translation-invariant (The interactions depend only 
on relative positions.), the potential energy will not change. 

The same statements can be carried over to quantum mechanics. We can 
relate the wave function of our system at rest to the wave function when it 
has velocity v by 



N 



exp 



(mi/h)^v- fi 

i=l 



where the index i ranges over the atoms of the system. Since the velocity 
enters only the phase factor, the expectation value of the potential energy 
will not change. 

The momentum and energy eigenvalues of the system in the two reference 
frames are: 

= P + Nmv 
^(^) = E + P ■v + Nmv^/2. 

Corresponding expressions for the ground state are: 

P^^ = Nmv 

= Eo + Nmv^/2. 

For the state in which a single quasiparticle (p, e) is present: 

= P + Nmv 
Ef = Eo + e{p)+p-v + Nmv''/2. 

Thus the change in momentum and energy associated with the single quasi- 
particle is 



17 



= p 

which is our equation ^ of section || 
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Figure captions 



Figure |T| Dispersion curve for elementary excitations in liquid helium, cal- 
culated from the polynomial fit in reference ||^. The roton region is 
near the minimum. 

Figure |2| Surfaces of constant energy in momentum space for e = A (dot- 
ted), and a higher energy (solid). 

Figure |^ Momentum space diagram showing initial and final momenta cor- 
responding to reflection of rotons by a force in the —y direction. 

Figure |^ Coordinate space trajectories corresponding to momenta in the 
three regions shown in the previous figure. Notice the different scales 
in the two directions. 

Figure |5| Energy of an elementary excitation, and an atom outside the liq- 
uid, and the "Refractive index" plotted against momentum/ h. 

Figure |^ Calculated roton trajectories near a quantized vortex line, located 
at the origin. The curves are calculated for an incident roton with 
momentum l.lpo incident from the right (solid lines), and an R~ roton 
with momentum 0.9po incident from the left (dashed lines). 

Figure |^ Evaporation of a He atom by an R~ roton. 
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